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Path integrals similar to those describing stiff polymers arise in the Helfrich model for membranes. 
We show how these types of path integrals can be evaluated and apply our results to study the ther- 
modynamics of a minority stripe phase in a bulk membrane. The fluctuation induced contribution 
to the line tension between the stripe and the bulk phase is computed, as well as the effective inter- 
action between the two phases in the tensionless case where the two phases have differing bending 
rigidities. 



I. INTRODUCTION 

Recently their has been much interest in the effective interactions between components of membranes of different 
composition. These effective interactions can be direct basic interactions such as electrostatic and van der Waals 
forces. However the fact that the membrane fluctuates also leads to effective interactions which are due to how the 
different components, or inclusions, alter the membrane fluctuations. Coarse grained models, based on the Helfrich 
model[l], 0, of multicomponent membrane describes the membrane in terms local mechani cal prop erties, for instance 
the bending rigidity Kb, the Gaussian rigidity K g or the spontaneous curvature [13, [0, [Hi 0, [13, [IB]- The 
effective, fluctuation mediated, interaction between regions of differing rigidity can be computed using a cumulant 
expansion giving the effective pair-wise component of the interaction between two regions. This term is of order 
^ K b/g an d is the analogue of the pair-wise component of van der Waals forces. However when | <5«b / 9 1 is large then 
this two-body, or dilute, approximation will break down and a full iV-body calculation is needed. We should expect 
the dilute approximation to break down reasonably frequently as experimentally measured values of Kb for commonly 
occurring lipid types vary from 3 to 30fcsT. 

In this paper we show how the full iV-body calculation for a system with spatially varying rigidity and elasticity, 
can be carried out for stripe geometries of the type shown in Figure [TJ Within this geometry we can compute 
the contribution to the line tension between the two phases due to membrane fluctuations. This can be seen as a 
renormalization of the line tension already present due to basic interactions such as van der Waals, electrostatic and 
steric forces. In addition, in some cases, we can evaluate the effective interaction between the two interfaces as a 
function of their separation I. 

They key point in our calculation is that we convert the usual functional integral into a path integral where the 
direction in which the physical parameters change is treated like a fictitious time variable within the path integral 
formalism. The authors have already applied this approach to electrostatic problems where it has proved to be efficient 
for carrying out computations for films pj] , interfaces [H and in cylindrical geometries such as lipid-tubules [H . 

The paper is organized as follows. In section [TT] we describe the model and show how it can be analyzed using path 
integrals which are mathematically identical to those arising for stiff polymers. In section Hill we present the results of 
our computations. We start with an analysis of bulk homogeneous membranes and show how some standard results 
can be recovered using our path integral formalism and then how the method can be applied to a striped membrane 
system. The formalism is then used to compute the membrane fluctuation-induced contribution to the line tension 
between two phases. In section HVl the fluctuation-induced Casimir force in a striped system is calculated for typical 
physical situations. In section|V]we conclude with a discussion of our results. The generalized Pauli-van Vleck formula 
used to evaluate the path integrals is described in detail in Appendix [XJ Some aspects of this approach have been 
previously described in the literature [1, H, Q, however our approach is slightly different and self-contained. Hence for 
the sake of completeness (and because the results seem to be relatively unknown) we include this detailed description 
of the approach. 
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FIG. 1: Schematic diagram of striped membrane configuration with mechanical parameters (rigidity and elasticity indicated. 



II. THE MODEL 



In the Monge gauge the Helfrich model for a membrane whose height fluctuations are denoted by h is given by 

H=\( d 2 x \k{V 2 K) 2 +ii{Vh) 2 ] . (1) 
2 Ja p 1 J 

In this model we have neglected the Gaussian rigidity n g and so for ease of notation we denote the bending rigidity 
simply by k. In this form of the Helfrich model there is no spontaneous curvature and it is implicitly assumed that 
the fluctuations of h are small. The term fi can be interpreted as a local elastic or surface energy. When /i is constant 
it can be interpreted as a surface tension. The integral in x = (z,y) in Equation (TT]) is over the projected are of 
the membrane A p . The physical area of the membrane is larger than A p due to fluctuations and we denote it by 
A = A p + AA, where AA is the excess area due to fluctuations. In the limit of small height fluctuations (i.e. to 
quadratic order in h) the excess area is given by 

AA = - f d 2 x {Vhf . (2) 

2 JA V 

The canonical partition function can be written as a functional integral over the height field h 

Z = J d[h}exp(-(3H), (3) 

where [3 — 1/ksT, T is the canonical temperature and fcs Boltzmann's constant. 

If the mechanical parameters n and /i only vary with the coordinate z then we can express h in terms of its Fourier 
decomposition in the direction y writing 

v ^ k 

We have imposed periodic boundary conditions in the y direction and thus have k = 2itn/L, where L is the width of 
the system and n is an integer. Note that we are assuming that the interfaces between differing phases are straight and 
and thus lie at constant values of z. This is a realistic assumption if the bare (in the absence of height fluctuations) 
line tension 70 is positive and large. We will in fact see that the renormalization of the line tension 7 between phases 
due to the height fluctuations is positive and thus this assumption remains valid (and is in fact reinforced) by height 
fluctuations. 

In the limit of large L the sum over modes can be written as 



k 

We now find that the Hamiltonian decomposes as 



k 



(6) 
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where 

it J d 2 h(z,-k)d^h(z,k) \ 2 ^( dh(z,-k)dh(z,k) \ ,, 2 - - 

H k = -ldzn{z)\ — l+(2K(z)fc +/i(z))l — — — I + (/i(z)fc + K(z))h(z, k)h{z, -k). 

(7) 

The field h is real and so we have the relation h(z, k) = h(z, —k). The full functional integral for the partition function 
Z can then be written as 



where 



z= JJe*. 



fe>0 



(8) 



9* = / d[X]expf~ jf 



dt 



a 2 (k,t) 



d_X\ 2 
~d7J 



i£) +<*i(M)(^) + a (M) a *t 2 



(9) 



In the above, the coefficients are given by 



a 2 (k,t) = /3 K (i) 

ai (fc,i) = (i{2n{t)k 2 + fi(t)) 

a (k,t) = /3(^(i)fc 2 + *(£)), 



(10) 



and Z is the length of the system (in the z direction). The above path integral is that arising in an elastic model of 
a semi- flexible polymer 0] (in one dimension) where k = a 2 /{3 is the rigidity of the polymer, the term proportional 
to a\ represents the elastic energy and the term proportional to ao represents an external harmonic potential. In this 
model the length of the polymer is not fixed, in contrast to the worm-like chain model where the magnitude of the 
tangent vector is fixed. The method of evaluation of the above type of path integral is given in Appendix [S] and we 
find that it takes the form 



K(X,Y;t) = (2tt)-- [det (B(t)))}* 

exp (-^X T A T (t)X - ^Y T A F (t)Y + X T B(t)Y 



(11) 



where the initial condition vector is X = (X, U) with X = X(0) and U — dX/ds\ s= o, and the final condition vector 
is Y = (y, V) with Y = X(t) and V — dX/ds\ s=t . When the coefficients a k are independent of t the classical action 
can be written as a combination of surface terms (using the equation of motion) as 



S cl (X,Y) = - \a 2 



dXd 2 X 



ds ds 2 



a 2 



X 



.d 3 X 



ds 3 



a i 



X- 



,dX 
ds 



(12) 



The above expression is in general rather complicated but can be used to determine the matrices Ap, Aj and B. Also, 
when the coefficients are independent of time, the time reversed trajectories have the same weight as the original 
trajectories. Thus the path integral going from (X, U) to (Y, V) in time t has the same value of the path integral 
going from (Y, — V) to (X, —U) in time t. Mathematically this means that 



A F = SArf, 



(13) 



and 



B 1 = SBS, 



(14) 



where 



S 



1 

-1 



(15) 



In the limit of large t the propagator K should factorize, as it is dominated by the lowest eigenfunction, and so we 
should find that B(t) — > as t — > oo. This can be verified by explicit calculation in the cases considered here. 
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III. CALCULATIONS FOR BULK AND STRIPED SYSTEMS 



To start with we will show how the path integral formalism introduced here reproduces some standard results 
concerning bulk systems. We consider a bulk system of projected length I and projected width L; we thus have a 
projected area A p = LI. Periodic boundary conditions are imposed in both directions z and y. The free energy is 
given by 

F = -k B Tj2M®k), (16) 

k 

where 

9 fc = J dXK k (X,X,l). (17) 

Using Equation (|A17j) and the notation developed in the appendix we find 

e fe = det (B(k, I)))* det (A/(fc, I) + A F (k, I) - 2B(k, . (18) 
The classical equation of motion in this case has solutions 

X(t) = acosh(pt) + 6sinh(p<) + ccosh(gi) + dsmh(qt), (19) 

where 

p = k, (20) 

and 

q = (k 2 +m 2 )?, (21) 

where we have defined 



2 M 

m = — , 



(22) 



and so m is an inverse length scale. The expressions for Ap Aj and B can be computed using computer algebra but 
they simplify in the (thermodynamic) limit I — * oo. We define 



A* I/F = km A I/F , (23) 

' l—>oo 



and find that 

.4 



*=p K ( pi(p+i) p<i y (24) 

\ pq (p + q) J 

We also find that for large I B(l) — * 0, and 

det (B(l)) « pq(p + qf exp (-(p + q)l) . (25) 
The extensive part of the bulk free energy is thus 

F k B T 



A p 2tt 



dk 



k+{k z + m z )i , (26) 



where we have introduced the ultra-violet cut-off length scale a which corresponds to the lipid size. The infra-red 
cut-off scale (where needed) is given by L (the lateral size of the system). The excess area of the system is given by 

dF 

AA=-, (27) 

and in the tensionless limit where fi = we find the well known result 

AA_k B T, (L 
A„ Attk 



In - . (28) 
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In the case where fi^Owe find that 



M = M: sinh -i(M. 

A„ Ann \amJ 



(29) 



This gives 



AA 



Ann V am 



(30) 



when a <C 1/m. 

For a striped geometry where the length of the bulk phase is Z , and large and that of the minority phase is I we 
find that for this composite striped (hence the superscript s in what follow) system we have, as Iq — * oo, 



ei s) (Uo) = J dXdYK k (X,Y,l)K { °\x,Y,l ) 

det(B< >(Z ))] 5 [det(B(l))]* [det (a<?>* + A,(0) 

(0)* 



det 



(4 0) * + a f (o - b t (o(4 0) * + aj(o) _1 5(o) 



(31) 



where the superscript (0) refers to the bulk phase and the absence of this superscript refers to the minority phase. In 
the limit I — > oo the above expression simplifies giving 



e { k s \l,lo) « [det(B(°)(/ ))l 2 [det(B(0)]' [det(4 0) *+^) 



(32) 



where we have used Equation (| 13|) . In order to compute the free energy cost of the interface between the two phases 
we subtract the separate bulk free energies for large systems of size lo, corresponding to the bulk, and of size Z, 
corresponding to the minority phase, from that of a large striped system composed of length Iq of the bulk phase and 
I of the minority phase. This free-energy difference is 



AF = -fc B T^ln 



ej B) (Uo) 
ei o) (z o )e fe (0, 



(33) 



which gives 



det(4 0) * + A$>*) det{A} + A* F ) 
det(4 0) * + A* F y 



|(0) 



(34) 



The above expression is in general quite complicated but when k and Ko are non-zero, then at large k the eigenvalues 
q (given by Equation (|21jl ) in the stripe phase becomes asymptotically equal to go, the corresponding eigenvalue in 
the bulk phase. We thus find that the sum in Equation ([B"4"]) is ultra-violet divergent and is dominated by the term 



AF=^ln 



1 - A 2 /4 



(35) 



where 



A = 



K — Kq 
K + Kq 



(36) 



We may interpret this result as the existence of a height fluctuation induced line tension 7^ f between the two phases 
(note the factor of a half as there are two interfaces) given by 



7hf 



2a 



In 



1 - A 2 /4 



(37) 



Thus the dominant contribution to the fluctuation induced line tension between the two phases comes from the miss- 
match in their bending rigidities. We also remark that it does not depend on m and is only dependent on k and kq 
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FIG. 2: The function 1(A) denned in Equation pg|l. where A = (« - Ko)/(k + k ). Note that 1(A) = /(-A). 

through A 2 , which is a symmetric function of the two rigidities. As an example the fluctuation induced line tension 
between two phases whose rigidity differs by a factor of 10 has an energy of about 0.5 fcsT per lipid at the interface. 

The correction terms to jhf are UV convergent and, after some manipulation, 7 can be expressed as a power series 
in ma/ir as 



Ihf 



1 /l-A 2 /4 

2 ln TTaT 



™ 1(A) -L 

7T v ' 128 



A 2 



(1- A 2 /4) 



1 /max 4 1 /ma 
3 \ V) ~ 5 V~ 



(38) 



J(A) is shown in Figure [2] for A > and is a non-negative function of A with 1(0) = and 1(A) = /(—A). From 
Figure [H we see that 1(A) has a maximum value I (00) ~ 0.04. For /z = 10 _2 -/V/m, k = 25/csT ~ 10~ 19 J and 
a = 10~ 9 to, we find ma — y/[fJ./ n)a ~ 0.32 and so ma/n ~ 0.1. The correction to the leading term due to non-zero 
/x is thus expected to be certainly less than 0(1%). 



IV. 



APPLICATIONS 



In this section we discuss the application of the theory developed above to two cases in the system with a stripe 
as shown in Figure [T] These cases are distinguished by the values of the masses in the two regions, too = ^/(/io/fto), 
to = y^l/i/n), which control the relationship of the surface to bending energies. In general, the boundary conditions 
satisfied by the system will vary and will determine the precise way in which our formalism is applied and the form 
taken by the relevant Helfrich action. We consider two cases that might be thought of as extreme situations and are 
chosen to show how the results are markedly different depending on the exact situation. We concentrate on computing 
the Casimir force across the stripe which can be interpreted as a force between the opposing interfaces. The Casimir 
free energy, Fc(l), is therefore normalized to F = in the limit I — > 00. We have 



F c (l) = F(Uo) - I™ F(IJ ) \ 



i+i =constant > 



(39) 



where 



F(Uo) 



fcsT^ln 



(©i s) (Mo; 



(40) 



and @]?\l,lo) is defined in Equation (|3"Tj) . It is understood that the total volume of the system is held fixed by 
imposing I + Iq = constant, and that Iq is large compared with any system-specific length scale. 
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FIG. 3: The function C(A) defined in Equation flU}, where A = (k - k )/(k + k ). Note that C(A) / C(-A). 



A. mo = m = 

This case corresponds to an untethered membrane which is tensionless, as is the case for a membrane in the presence 
of the various lipid species in solution. Upon a change in the physical area of the membrane, A, lipid molecules can 
leave or enter meaning that any area change costs no free energy. 

Then we have /io = /i = 0. In this case, the choice for the general solution to the classical equations of motion is 
not given by Equation (fT9|) but by 

(41) 



X(t) = acosh(pt) + bsinh(pt) + ctcosh(pt) + dts'mh(pt) 
The method follows the manipulations of section IIII1 and appendix [A] We find that 

k B T 



In (1 + a 2 (l, A)e' 2kl + a 4 {l, A)e 



4k-l\ 



(42) 



where 



020, A) 
a.4(l, A) 



(1- A 2 /4) 

9A 4 
16(1- A 2 /4) 2 



2;2 f l-A/2 
1 + A/2 



(43) 



We note that all terms are invariant under A — > —A except the first term in 0,2(1), and hence the Casimir free 
energy is not invariant under this transformation in this case. On dimensional grounds we have 



C(A)k B T 
I 



and C(A) is shown in Figure [31 Since C(A) < for A / 0, the Casimir force is attractive and is given by 

d C(A)k B T _ C(A)k B T 



/c(J,A) = - 



dl 



r- 



(44) 



(45) 



with |C(A)| £ 0.4. 
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FIG. 4: The function G(T) defined in Equation gHJ, where T = (fj, — (J,o)/(fi + Mo)- Note that G(T) = G(-T). 



m , m > 



This case is for non-zero /iq and /i. The projected area is constant but any change in the area results in a free 
energy change. In practice, the typical observed values of \i give ma ~ 0.3 (see section IIII[) and consequently the 
width of the stripe, I, satisfies ( < As shown in section UTT1 the free energy has a contribution that corresponds 

to the energy of the interfaces between the stripe and the bulk medium and includes the UV divergent part of the 
sum over modes. The remaining, ^-dependent, terms are UV convergent and are cut-off but exponential factors on a 
scale k ;$ Thus, the Casimir free energy excluding the interface energy gets contributions only from low mode 
number k <C 1 /m and so we can approximate the eigenvalues in Equations (|20p and (|2ip by 



P 



m 



(46) 



Also, only leading terms in e ql will survive; the others will be suppressed by factors of e~ ml . Following the derivation 
of the previous subsection, we find a result that is independent of Kq and k, as one might expect on dimensional 
grounds for fj, » njl 2 . We get 



Fg(IX) = ^£i n (i_r 2 — " 



■') 



(47) 



where T = (/x — ^o)/(M + Mo)- The leading corrections are suppressed by the factor (1/ml) 2 . For the values of m 
considered this is of order (a/l) 2 ~ 1/A^ 2 , where N is the number of lipid molecules across the strip. Similarly to 
before we may write 



F G (l,T) 



G(T)k B T , / 1 



O 



(48) 



G(r) is symmetric under T — > — T, and is shown in Figure^ We see that |G(r)| ^ 0.2 and, being negative for V ^ 0, 
gives rise to an attractive Casimir force of maximum magnitude 0.2kgT /I 2 . 

The more general case where m ~ 1// is very much more complicated and the expressions cannot be presented here 
but need to be investigated computationally in the three parameter fi, k, I space. We have verified that this analysis 
is feasible but complicated and so we have chosen not to present it in this paper. However, any more general result 
for the Casimir force will interpolate between the two extremes presented here and so we would expect a leading 
contribution to behave like —cksT/l 2 with c J§ 1. 
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V. CONCLUSION 

In this paper we have shown how the formalism developed in earlier work 0] can be applied to the more general case 
of higher derivative Gaussian energy functions such as apply to the path integral analysis of stiff polymers and the 
Helfrich model for membranes. Some aspects of this approach have been previously described in the literature [H, [H, [|| . 
However, our approach is slightly different and self-contained. Hence, for the sake of completeness (and because the 
results seem to be relatively unknown) we have included a detailed description of this approach. In particular, we 
have shown how to generalize the Pauli-van Vleck formula for the evolution kernel of all theories of this type. 

As a model system we have considered a toroidal lipid membrane with one very large circumference and the other 
finite of length L, with a stripe of width / wrapped around the finite circumference and of different, minority, lipid 
type to the bulk, majority, type. This is shown schematically in Figure [TJ This geometry imposes periodic boundary 
conditions on the system. We have shown how to compute, in general, two important energies in this system, namely 
the energy, or line tension, associated with the lipid-lipid interface and the Casimir force between the interfaces, as a 
function of width I. An important controlling parameter, to, has the dimensions of a mass and is given by m 2 = fi/ k. 
We have presented these calculations explicitly for the cases where to = and ma ~ 1, I > where a is the 

inter- lipid spacing. These correspond, respectively, to the cases where the actual area, A, or the projected area, A p , 
is conserved. In the latter case, we calculate the mean excess area of the system AA/A p in section ITTT1 In both cases, 
the interface energy is positive and the Casimir force attractive as can be seen from Equations (|35|) , (|44p . (|48|) and the 
associated figures. Our general result is that in appropriate dimensionful units the energy coefficient is cksT where, at 
maximum, c ~ 1. This is to be compared with the natural bending rigidity with lies in the range 5ksT < n < lOOksT. 

The general case where I ~ 1/to is complicated and long, and although we have the results we have not presented 
an explicit analysis in the /i, k, I parameter space because of the complexity. However, there is no computational or 
algebraic impediment to carrying this out. 

In terms of relevance to the physical system we might consider two scenarios in the two-lipid model discussed here. 
Either the minority lipid can be dissolved in the majority lipid to form a homogeneous phase for the mixture, or the 
minority lipid can precipitate out of solution and form a pure minority phase, the stripe in our idealized case, within 
the pure majority phase. Which situation is stable is, of course, decided by a competition between the free energies 
of the configurations which is in turn dependent on the boundary conditions imposed. However, it is clear that the 
attractive Casimir force will tend to reduce the stripe width I, presumably by evaporation of minority lipid from the 
interface into solution. The interface energy is constant throughout such a process but will always tend to minimize 
the interface length. A stability analysis, however, requires a computation of the free energy of the mixed phase which 
our calculation does not address. However, as has been discussed in (l6j . the suppression of lipid mode fluctuations by 
confining the membrane in a stack will change the free-energy of both configurations and so can affect their stability; 
an effect which can be analyzed by our methods. 



APPENDIX A: THE GENERALIZED PAULI-VAN VLECK FORMULA 

In this appendix we show how generalized quadratic path integrals can be evaluated giving a generalization of 
the Pauli-van Vleck formula. The treatment is very close to that of @ and is based on the Chapman-Kolmogorov 
decomposition of the path integral. We consider the following path integral 

tf(X,Y;t)= / X(t) Y d[X] exp(-5[A]), (Al) 
Jx(o)=x 

where S is a quadratic action which will have the general form 

J{> k=o v y 

In general the coefficients at can be time dependent but for the problems related to membranes studied here we will 
only require the results for ak constant. The usual Wiener measure occurring in path integrals has N = 1 and the 
corresponding path integral is that for standard Brownian motion, or a free particle, with clq = 0. If ao m this 
case, then the path integral corresponds to that of a simple harmonic oscillator with ao = tow 2 and a\ — m thus 
relating the coefficients i = 0, 1 to the mass to and frequency of the oscillator. The path integrals arising in 
section [XT] are, of course, for the case N = 2 which, as mentioned above, also arises for the path integrals of stiff or 
semi-flexible polymers. We also refer the reader to the approach of [3| which is based on an eigenfunction expansion 
method for the case N = 2. 
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Now one must state how the initial and final points of the path integral should be specified. The presence of 
the term (d N X/ds N ) means that the paths that contribute to the path integral are ones where the derivatives 
d N ~ 1 X/ds N ~ 1 and lower must be continuous. The path integral should therefore be specified in terms of the vector 
X = [X,X^ X \X^ ■ ■ ■ X^^ 1 ^) where we have used the notation X^ = d k X/ds k . We can now decompose the path 
integral using the Chapman-Kolmogorov formula 

K(X, Z;t + t') = J dYK(X, Y; t)K(Y, Z; t'). (A3) 

This decomposition ensures the continuity of the path X(t) up to its TV — 1th derivative. 
The classical path is given by the one that minimizes the action: 

™ 0, (A4) 



SX(s) 

with the boundary conditions on the end points X(0) = X and X(t) = Y. This gives a total of 2N boundary 
conditions (N from each end). The equation for the classical path can be written as 

ds 'sxM^) Xcl(s ' ) = ' (A5) 

which is a linear differential equation of order 2N. For instance, when the are constant it reads 

N rl 2k 

E(- 1 ) fc °*^fc^( fl ) = - ( A6 ) 

fc=0 

There are thus 2N linearly independent solutions to this equation and their coefficients are linearly related to the 2N 
conditions for the end points. The classical action is a quadratic form in the initial and final condition vectors X and 
Y and we can write 

S cl (X, Y) = i [X T Aj (t)X + Y T A F (i)Y - 2X T B(t) Y] , (A7) 

where we have used the subscripts I and F to denote the initial and final coordinates. We now write the path X(s) 
as X(s) — X c i(s) + x(s), where the boundary conditions imply that x(0) = and x(t) = 0. The path integral can 
now be written as 



K(X,Y;t) = exp(-S c; (X,Y))^ d[x]exp(~J dsds' x(s' ) . . 



S 2 S 
6X(s)SX(s'Y 

Q(t)exp(-S d (X,Y)), (A8) 



where we formally write can write 

Q ^ = det { mBm ) 2 ' °- s ' s '-" (A9) 

The above functional determinant can be evaluated using an eigenfunction expansion, however for higher order 
operators this quickly becomes impractical. Instead, we return to the Chapman Kolmogorov formula Equation (|A3D 
and pursue its consequences using the formal result Equation (|A8[) . Explicitly carrying out the intermediate integration 
over Z, we find that 

K(X,Z;t + t') = (2TT)^Q(t)Q(t r )det(A I (t)+A F (t')y i x 

cxp (~X T [Ai(t) - B(t)(Ai(t') + A F (t))- 1 B T (t)] X 

exp (~Z T [A F {t>) - B T {t'){A I {t') + Apit^Bit')} Z 
cxp (X T B(t)(A I {t') + A F (i))- 1 J B(t')Z) , (A10) 



Now comparing the quadratic forms and prefactors we find the following relations: 

A I (t + t>) = A I {t)-B{t){A I {t')+A F {t))- 1 B T {t) 
A F (t + t') = A F {t')-B T {t l ){A I {t l ) + A F (t))- 1 B{t') 
B(t + t') = B(i)(A / (t') + A Ji .(i))- 1 .B(t') 

Q(t + t') = (2tt)* Q(t)Q(t') det (A/(t) + A F (t'))~ h . 



(All) 
(A12) 
(A13) 

(A14) 



In Q it is pointed out that relation Equation (|Al4j) above can be used to derive a differential equation for Q. However, 
a more rapid way of finding Q is to note that taking the determinant of both sides of Equation (j A13|) gives 

det (B(t + t')) = det (B(t)) det {B{t')) det (Ai(t) + A F (t')) _1 , (A15) 
and using this relation we find that by direct substitution into Equation (|A14p that the solution for Q is 

Q(t) = (2w)-% [det(B(t))]*. (A16) 

The generalized form of the Pauli-van Vleck formula may thus be written in the familiar form (for N = 1 and its 
generalization to higher dimensions) 



K(X,Y;t) = (27r)--det 



dS cl 



exp (S cl (X,Y)) 



(A17) 
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